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Рассматривается горение твердых ракетных топлив и порохов, разлагающихся на газ посредством пиролиза. Показана возмож- 
ность гашения горения гармонически меняющимся давлением. Причиной тому является то, что при определенных условиях ма- 
лые колебания давления, как правило, вызывают большие по амплитуде изменения скорости горения. Погасание происходит 
при снижении температуры поверхности топлива ниже критической величины. Проведено исследование для случаев экзотер- 
мической и эндотермической реакции пиролиза. 


Введение 

Практика разработки и эксплуатации ракетных 
двигателей на твердом топливе встречается с рядом 
разновидностей нестационарного горения, проис- 
ходящего при переменном давлении. Изменение да- 
вления может быть вызвано штатными и нештатны- 
ми ситуациями. К нештатной ситуации относится 
неустойчивость горения, вызванная конструкцион- 
ными особенностями двигателя, внутренними фи- 
зико-химическими процессами, обеспечивающими 
горение топливного заряда, внешним воздействием. 
Такой вид нестабильного горения в ракетных двига- 
телях подразделяют на низкочастотную и высокоча- 
стотную неустойчивость [1, 2. С. 104-129]. Вторую 
из них называют еще акустической неустойчиво- 
стью или, реже, резонансным горением. 


Низкочастотная неустойчивость. В процессе 
горения топлива наблюдалось, что когда давление 
в ракетном двигателе падает ниже какой-то крити- 
ческой величины, давление в камере сгорания мо- 
жет внезапно упасть до атмосферного, и горение 
заряда прекращается. Иногда, однако, спустя неко- 
торое время, от долей секунды до нескольких се- 
кунд, вновь происходит воспламенение, и наступа- 
ет новый период нормального горения. Этот цикл с 
частотой -1...10 Гц может многократно повторять- 
ся, приводя к ряду вспышек. Высокочастотная 
(і акустическая ) неустойчивость наиболее часто на- 
блюдалась в трубчатых зарядах, хотя она встречает- 
ся и в случае зарядов другой конфигурации. Горе- 
ние при такой неустойчивости характеризуется 
внезапным ростом давления в ракетных двигателях 
и колебаниями с частотой порядка 10 3 ... ІО 4 Гц. 
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К штатным изменениям режимов горения от- 
носится запланированный сброс (или подъем) да- 
вления в камере сгорания с целью перевода двига- 
теля на режим с другой тягой. Смена скорости го- 
рения в двигателе может быть осуществлена и дру- 
гими способами: комбинированием различных ви- 
дов топлива, имеющих разные скорости горения 
при одних и тех же внешних условиях, изменением 
критического сечения сопла и др. 

Кроме этого, при создании реактивных старто- 
вых двигателей, предназначенных для выброса 
торпеды или ракеты из пусковой установки, стоит 
проблема обеспечения заданной скорости прихода 
газа путем изготовления топливных зарядов спе- 
циальной геометрической конфигурации [3]. Но и 
здесь неизбежно происходит изменение давления в 
камере в процессе работы двигателя. 

Для разработчиков ракетных двигателей при 
указанных выше факторах необходимо знать зако- 
номерности горения топливного заряда. Наиболь- 
ший интерес представляет нестационарная ско- 
рость горения и как функция давления р и времени 
і. Эту зависимость и(р, 1) можно найти разработкой 
соответствующей физико-математической модели 
горения топливного заряда и последующим ее ана- 
лизом. Для предупреждения возникновения высо- 
кочастотной неустойчивости необходимо еще под- 
держивать заданные пределы изменения акустиче- 
ской проводимости С, топлива. 

Настоящая работа посвящена расчету неста- 
ционарной скорости горения и акустической про- 
водимости топлива, установлению закономерно- 
стей их поведения при различных формах зависи- 
мости давления от времени. 


Математическая постановка задачи 


В одном из простейших вариантов модель горе- 
ния твердого ракетного топлива или взрывчатого 
вещества с его разложением по механизму пироли- 
за математически может быть представлена следу- 
ющим образом [4, 5]: 
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О начальных условиях для ур. (1) будет сказано 
ниже. Здесь Я, - коэффициент теплопроводности 
материала топлива; Т с — температура в твердой фа- 
зе, Т с (х=х ! )=Т ! - температура поверхности разложе- 
ния; У, Б - массовая концентрация (доля) и коэф- 
фициент диффузии реагирующего вещества; Я 0 - 
предэкспоненциальный множитель в законе Арре- 
ниуса; Е - эффективная энергия активации хими- 
ческой реакции в газовой фазе; Т - температура га- 
за; К - универсальная газовая постоянная; Я, - га- 
зовая постоянная; Я - коэффициент теплопровод- 
ности газа; 0 - суммарный тепловой эффект хими- 
ческой реакции в газе; у— показатель адиабаты; I - 
тепловой эффект разложения топлива на газообраз- 
ные компоненты, положительный знак которого 
указывает на эндотермическую реакцию, отрица- 
тельный знак - на экзотермическую. 

Область х<х, занята топливом. На поверхности 
х=х 5 происходит его разложение по механизму пи- 
ролиза на газовые компоненты, состоящие преи- 
мущественно из горючего и окислителя. В даль- 
нейшем эта смесь химически реагирует с тепловым 
эффектом, составляющим основной запас химиче- 
ской энергии топлива или взрывчатого вещества. 

Ур. (1) справедливы для установившегося режи- 
ма горения, поэтому в них отсутствуют начальные 
условия. Массовая скорость распада топлива от, 
описывается законом Аррениуса т~р ѵ,і ехр(-Е с /2КТ) 
с энергией активации Е; ѵ 0 =соп§1. 

При формулировке математической модели те- 
плоемкости твердой и газовой фаз полагались оди- 
наковыми: с=с } . В дальнейшем изменением моле- 
кулярного веса газа пренебрегаем: і? 8 =сошІ, 
с,,=сошІ. Также зависимости коэффициентов диф- 
фузии и теплопроводности от температуры и да- 
вления представим для удобства как [6] 

0~Т 2 /р, Я ~ Т. (2) 


Выражения (2) хорошо согласуются с теорети- 
ческими и экспериментальными данными. Прис- 
ваивая верхний индекс 0 к стационарным значе- 
ниям символов, переходим в системе ур. (1) к без- 
размерным величинам и лагранжевой координате д, 
согласно выражениям: 
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Переход к лагранжевой координате позволяет 
отделить гидродинамическую часть задачи от диф- 
фузионно-тепловой части. В результате несложных 
преобразований приходим к следующей формули- 
ровке задачи: 
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Задача (3) упрощается, если /3— >0, так, что 

| Ж 0 аф(-^^ = С0П5І. 


Тогда выражения с экспонентами, описываю- 
щие химическую реакцию в газе, можно перенести 
в граничные условия, введя поверхность горения у. 
На этой поверхности массовая скорость горения 
т 2 ~ехр(-Е/2КТ). 


В такой упрощенной постановке данная задача 
рассматривалась в работе [5]. Ее стационарное ре- 
шение дается формулами 

<^<0:0° =0 о + (1-0 о )ехр<§; (4) 

О<^<<^:0°=0 о -/ + (1-0 о + /)ехр^ 

Е = 1 - (1 - Т°) ехр(^ / Ье); 

I > : 0° = 0 О + ^ - /, Т = 0; 5° = 1, 


где У° - массовая концентрация горючего компо- 
нента вблизи поверхности разложения топлива; 
0 4 ° - безразмерная температура пламени: Ѳ^=Т^/ТУ 


Результаты исследования и их анализ 

Задача (3) решалась численно по неявной схеме 
со вторым порядком точности по шагу простран- 
ственной переменной и первым порядком точно- 
сти по временному шагу. В качестве начальных 
условий были приняты решения (4). 

Константу ѵ 0 полагаем пока равной нулю. Это 
означает зависимости скорости разложения топли- 
ва непосредственно только лишь от температуры 
поверхности топлива. Существуют эксперимен- 
тальные данные [7], хотя и не совсем надежные, но 
все же указывающие на отсутствие прямой зависи- 
мости скорости разложения от давления. Результа- 
ты моделирования нестационарного горения каче- 
ственно зависит от знака теплового эффекта Г. По- 
этому следует отдельно разобрать случаи экзотер- 
мической и эндотермической реакции на поверх- 
ности топлива. В настоящей работе рассмотрение 
ограничено низкочастотными колебаниями давле- 
ния. Именно в этой области наблюдаются наибо- 
лее интересные критические явления. 

Экзотермическая реакция. Анализ модели горения 
проводился на основе двух форм зависимости давле- 
ния от времени: синусоидальной и экспоненциаль- 
ной. Если состояние горения такое, что оно находит- 
ся вблизи границы устойчивого горения, то неболь- 
шие колебания давления могут привести к погаса- 
нию. Причем в случае синусоидальной зависимости 
давления от времени ? 7 = 1 +А, ) - 8 Іп(й>г) для каждой ча- 
стоты со существует свое критическое значение ам- 
плитуды А ц колебания, зависящее от остальных фи- 
зико-химических параметров. Аналогичная ситуация 
имеет место для экспоненциального сброса давления 
?7=1-Д?7+Д?7'ехр(-т/ г,,). Он характеризуется величи- 
ной перепада давления Д ?7 и временем г,, в течение 
которого величина давления уменьшается в е раз. От- 
метим, что в приведенных формах закона изменения 
давления для простоты записи опущено время за- 
держки действия нестационарной части. На приве- 
денных ниже рисунках сначала устанавливается ста- 
ционарное горение в течение времени задержки. 
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Мера удаленности состояния горения от грани- ный характер колебания (рис. 1, а) ярко выражен 

цы устойчивости определялась по величине фено- (В противном случае все кривые были бы синус - 

менологических коэффициентов г и к. Изначально оидальными). К тому же А ц близко к критическому 

их предполагалось находить экспериментальным значению 0,51. При более высоких амплитудах 

путем [7]. Но можно определить их теоретические происходит погасание (рис. 1, б). Его наступление 

значения, которые будут зависеть от принятой мо- характеризуется внезапным и быстрым падением 

дели горения. В нашем случае справедливы форму- скорости горения и температуры поверхности 

лы [ 8 ] (штриховая линия). Значения безразмерных пара- 

, 0 \2 метров следующие (вариант 1): 0 о =О,6; /=-0,23; 

г - Р' о - Е % ^ д=2,41; /3=0,06; /3=0,02; ^=3,8-10 8 ; /=0,14; к= 0,95. 

2-Р(Ѳ°) 2 Е с ут° ) ’ Исходные константы 0 О , /, < 7 , Д, /3, 5Г 0 подбирались 

о т такими, чтобы коэффициенты к, г были в пределах 

к = 0 ^ д ( 5 ) экспериментально установленных значений [7, 9]. 

2/3(0 ь )‘ 2К(Т Ь ) На рис. 1, 2 безразмерное давление становится от- 

Расчеты проводились для интервала частот личны м от единицы в моменты времени г>3. 

0,1. ..18 в единицах (и°) 2 /к с . Более высокие частоты, С увеличением частоты критическая амплитуда, 

как выяснилось, рассматривать нецелесообразно. как правило, тоже растет, но не значительно (рис. 2 , 

Во всяком случае, при о>>18 каких-либо новых ю=3). Здесь критическая амплитуда равна 0,442. 

критических эффектов не было обнаружено. Обращает на себя внимание постепенное зарожде- 

Пример поведения скорости горения В, темпе- ™ е нов °й моды колебания. На рис. 1 заметна свое- 

ратуры поверхности Ѳ, разложения топлива и кон- образная неглубокая впадина при 5*1,1. Эта впади- 

центрации У, горючего на этой поверхности приве- на П Р И частоте ю=3 (рис. 2) практически сформиро- 

ден на рис. 1 (для всех рисунков ст=0,01; Бе=1). вал „ ась в колебательную моду. Если период нели- 

Здесь ®=1. Ввиду большой амплитуды нелиней- неиного колебания на рис. 1 был равен примерно 



б 

Рис. 1 . Нелинейные колебания в системе: а) А^=0,4 вблизи критического значения амплитуды давления и наступление погаса- 
ния; 6) А, =0,52 
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Рис. 2. Нелинейные колебания в системе: а) А,=0,43 и критический переход ; б) А ц =0,443 к погасанию 
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т 


Рис. 3. Переходный процесс при амплитуде А п =0,41 и частоте 

4,9, то на рис. 2 он составляет 2, 6. ..4, 2. На профиле 
концентрации У, с увеличением частоты начинают 
появляться острые выступы (рис. 2, б). На стадии 
погасания температура снижается волнообразно, 
но неуклонно падает. При частоте й)=10 критиче- 
ская амплитуда равна 0,45. Основные качественные 
свойства переходного режима горения здесь такие 
же, как и на частоте ю= 3. Но на частоте ы= 1 5 
(рис. 3) погасание происходит не гладкими паде- 



3 4 5 т 6 


ю=15 

ниями основных динамических параметров, оно 
сопровождается сильными скачками скорости го- 
рения б, температуры Ѳ ( и концентрации У,. 

Здесь снова происходит снижение значения кри- 
тической амплитуды. Погасание происходит после 
того, как скорость горения или температура поверх- 
ности испытывает слишком резкий подъем и такое 
же падение. Причем на стадии падения появляется 
излом, после которого эти параметры испытывают 
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Рис. 4. Переходы стационарного горения к погасанию при превышении критических значений амплитуды давления: 
а) А =0,05, (0=3; 6) А =0,043, со=7; в) А =0,2; со=12 


еще большее снижение по величине. Если в пред- 
ыдущих случаях погасание произошло примерно за 
3 ед. времени, то теперь это время равно 1,5. 

Следующая серия расчетов (рис. 4) проведена 
для других исходных параметров горения (вари- 
ант 2): 0 о =О,3; /=0,23; <?=2,2; /1=0,08; /3= 0,04; 


А' 0 =3,3- ІО 5 ; г=0,13; Іс= 1,17. Действие нестационар- 
ного давления начинается в момент времени т=5. 

Если в 1-м варианте критическая амплитуда ме- 
няется слабо, то здесь заметен ее рост. Время погаса- 
ния, которое определяется грубо по характеру спада 
скорости горения, соответственно равно 3,5; 2,5; 1,3. 
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Технические науки 




а б 

Рис. 5. Погасание горения при повышении давления, Ар=-0,1; а) г, =7; 6) т ц =2 




Рис. 6. Погасание горения при эндотермической реакции: \г)=~0,22 ; т г =1/3 


Проведенное исследование для варианта 2 с эк- 
споненциальным изменением давления показало 
следующее: если перепад давления ниже значения 
Д?7«0,85, то система переходит асимптотически к 
новому стационарному режиму горения с повышен- 
ной скоростью горения В> 1. Если же выше 0,85, то 
скорость горения становится ниже первоначального 
значения 1. Также скорость горения понижается при 
повышении давления (т. е. при Д? 7 < 0 ). При этом от- 
носительно небольшое изменение давления приво- 
дит к погасанию горения (рис. 5). 

Поведение скорости горения в зависимости от 
характерного времени г, изменения давления гово- 
рит о слабой устойчивости рассматриваемого режи- 
ма горения. Это подтверждается тем, что при 
Д 77 = — 0,05 скорость горения, неуклонно понижаясь 
в среднем по времени, испытывает колебания с на- 
растающей амплитудой. В конечном счете, при 
прохождении через критическое минимальное зна- 
чение «0,3 скорость горения быстро спадает к нулю. 

Эндотермическая реакция. Такой тип реакции 
разложения топлива на газ рассматривался при 
следующих значениях параметров: 0 о =О,43; 
/=-0,23; ст=0,01; Ье=1; 4 = 2 , 6 ; Д=0,2; /3= 0,02; 
^=9,38- 10 9 ; 7=0,64; к= 1,82. Найденный критиче- 
ский перепад давления равен -0,215 (рис. 6 ). 


Если же давление меняется по гармоническому 
закону, то критическая амплитуда колебания ока- 
зывается равной критическому перепаду давления 
при экспоненциальном его изменении. Т. е. при 
протекании экзотермической реакции разложения 
на поверхности топлива для гашения горения тре- 
буется большая амплитуда давления. 

Выводы 

Анализ нестационарных режимов горения твер- 
дых ракетных топлив и взрывчатых веществ при 
переменном давлении показал, что: 

1 . Стационарное устойчивое горение может быть по- 
гашено при гармонически меняющемся давлении. 

2. Прекращение горения происходит при сниже- 
нии температуры поверхности топлива ниже 
критического значения. 

3. Имеется тенденция повышения критического 
значения амплитуды колебания с ростом его ча- 
стоты. 

4. При эндотермической реакции разложении то- 
плива отмечается повышенная устойчивость го- 
рения к колебаниям давления. 
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